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Comparing experimental data for high temperature cuprate 
superconductors with numerical results for electronic models, 
it is becoming apparent that a hopping along the plaquette di- 
agonals has to be included to obtain a quantitative agreement. 
According to recent estimations the value of the diagonal hop- 
ping t' appears to be material dependent. However, the val- 
ues for t' discussed in the literature were obtained comparing 
theoretical results in the weak coupling limit with experimen- 
tal photoemission data and band structure calculations. The 
goal of this paper is to study how t' gets renormalized as the 
interaction between electrons, U , increases. For this purpose, 
the effect of adding a bare diagonal hopping t' to the fully 
interacting two dimensional Hubbard model Hamiltonian is 
investigated using numerical techniques. Positive and nega- 
tive values of t' are analyzed. Spin-spin correlations, n(k), 
(n) vs yL, and local magnetic moments are studied for values 
of U /t ranging from to 6, and as a function of the elec- 
tronic density. The influence of the diagonal hopping in the 
spectral function A{k,uj) is also discussed, and the changes 
in the gap present in the density of states at half-filling are 
studied. We introduce a new criterion to determine proba- 
ble locations of Fermi surfaces at zero temperature from n(k) 
data obtained at finite temperature. It appears that hole 
pockets at k = (7r/2,7r/2) may be induced for negative t' 
while a positive t' produces similar features at k = (tt, 0) and 
(0, tt). Comparisons with the standard 2D Hubbard [t' = 0) 
model indicate that a negative t' hopping amplitude appears 
to be dynamically generated. In general, we conclude that it 
is very dangerous to extract a bare parameter of the Hamilto- 
nian it') from PES data where renormalized parameters play 
the important role. 



I. INTRODUCTION 

As experimental measurements of the properties of 
high temperature superconducting materials become 
more accurate, it has been observed that different cuprate 
compounds present slightly different properties. Among 
these qualitative differences are the presence of incom- 
mensurate magnetic correlations, the symmetry of the 
pairing state, shape of the Fermi surface (FS), as well as 
the behavior of the resistivity with temperature. 

Many models proposed to describe the cuprate super- 
conductors are particle-hole symmetric and this makes 



it impossible to distinguish between electron and hole 
doped materials. The addition of a diagonal hopping 
term to the one band Hubbard model has often been 
suggested as a way to handle the different properties 
between electron and hole-doped materials since such a 
term breaks particle-hole symmetry. Even restricted to 
the family of hole doped materials the inclusion of a di- 
agonal hopping has also been proposed in order to repro- 
duce experimental photoemission (PES) data and band 
structure calculations. According to these estimations, 
the diagonal hopping parameter t' seems material depen- 
dent and some of the suggested values are t' = 0.20 for 
Nd2CeCuOi, t' = -0.20 for LaiSrCuOi, and -0.45 for 

Recently, the effect of this t' hopping term has been an- 
alyzed in the weak coupling limit for the Hubbard model 
as well as for the t-J model. [^,|| It was reported 
that the Fermi surface of the U — t — t' model in the 
non-interacting case is in agreement with experimental 
measurements for appropriate values of t' , while the FS 
of the i' = model is not. The same occurs for the 
magnetic susceptibility and the sign change of the Hall 
coefficient with doping. Q However, it is not clear how 
a moderate or strong Coulombic interaction will modify 
the properties of weakly interacting electrons. In other 
words, it is not correct to deduce the value of a hare pa- 
rameter in a Hamiltonian using data from experiments 
where the renormalized parameter plays the key role. It 
is well known that strong correlations are fundamental in 
high Tc models, thus it is likely that the bare t' is quite 
different from the renormalized one. 

In this paper, we investigate the U — t — t' model in 
the intermediate coupling regime using quantum Monte 
Carlo techniques. We will compare the results with mean 
field calculations, as well as with the non-interacting limit 
and with experimental results. The paper is organized 
as follows; in Section II the model is introduced and its 
symmetry properties are discussed. Spin-spin correlation 
functions and incommensurability are studied in Section 
III. Section IV is devoted to the behavior of the density 
(n) as a function of chemical potential while local mag- 
netic moments are discussed in Section V. Angle resolved 
photoemission (ARPES) experiments and the shape of 
the Fermi surface are studied in Section VI. Dynamical 
properties appear in Section VII and the conclusions are 
presented in Section VIII. A mean field analysis is dis- 
cussed in the Appendix. 
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II. THE MODEL 

The U — t — t' Hamiltonian is given by 
H=-t ^ (4.Cj,, + h.c.)-t' (cUc„,, + h.c.)+ 

<ij>,(T <in>,(T 

U^(niT - l/2)(nu - 1/2) +^^11;,, (1) 

i i,(T 

where Cj^ creates an electron at site i with spin projec- 
tion (7, nio- is the number operator, the sum (ij) runs over 
pairs of nearest neighbor lattice sites, and the sum (in) 
runs over pairs of lattice sites along the plaquette diago- 
nals. U is the on site Coulombic repulsion, t the nearest 
neighbor hopping amplitude, t' the diagonal hopping am- 
plitude, and /i is the chemical potential. 

When t' = the above Hamiltonian is particle-hole 
symmetric, i.e., it remains invariant after the particle- 
hole transformation operation 

^ (-l)'ci.<.- (2) 

lit' > Q the above transformation maps the Hamiltonian 
in its particle- hole symmetric one but with < 0, and 
vice versa. This means that the results obtained for neg- 
ative (positive) t' below half-filling can be obtained from 
the results for positive (negative) t' above half-filling. In 
fact, this is the reason why to mimic the effect of electron 
and hole doping in the t-J model, (which cannot be doped 
above half-filling) both signs of t' are studied below half- 
filling. Then, it is clear that results above and below 
(n) = 1 with different signs for t' are redundant and in 
this paper we will study the Hamiltonian (Eq.(l)) only 
below half-filling with positive and negative values for t' . 
We will set t = 1 and |i'|=0.2 unless stated otherwise. 

III. MAGNETIC CORRELATIONS AND 
INCOMMENSURABILITY 

A spin density wave (SDW) mean field analysis (see 
Appendix) of the U — t — t' model predicts the existence 
of antiferromagnetism at half-filling for J7 > /7c- Q At 
\t'\ — 0.2 it can be shown that Uc = 2.1. Thus, within the 
SDW mean field, adding a t' term appears to decrease the 
strength of the spin-spin correlations. This result is not 
unexpected if we consider the large U limit. In this case 
a Heisenberg spin-spin interaction with J = At^ /U ap- 
pears between nearest neighbors sites and & J' — At''^ /U 
is generated along the diagonals. The J' interaction in- 
troduces frustration and tends to decrease the tendency 
to antiferromagnetism. In Fig.l-a we present quantum 
Monte Carlo (QMC) results for the spin-spin correlations 
C(r) = (S?S?_^j.}(— l)l''l as a function of distance, at half- 
fiUing, for U/t = 4, = and -0.2 on an 8 x 8 lattice 



at (3t = 6. Clearly, in both cases there is long range 
order (or, better, a correlation larger than the system 
size) but the strength of the antiferromagnetic correla- 
tions decreases when \t'\ increases. This is true for either 
sign of t' since due to the symmetry of the Hamiltonian 
at (n) = 1 the spin-spin correlations must be the same 
for both signs oit' . A similar effect is observed in Fig.l-b 
as U/t is increased to 6 at /3i = 4. 

Previous quantum Monte Carlo calculations on the 2D 
Hubbard model have indicated that short range incom- 
mensurate magnetic correlations develop as the system is 
doped away from half-filling. ||^,^ After these results were 
reported, neutron scattering experiments revealed short 
range magnetic incommensuration in La2-xSrxCuOi 
(LSCO). |£| The qualitative agreement between theory 
and experiment was remarkable. However, at a fixed dop- 
ing, the experimental peaks in the structure factor S(k) 
appeared shifted away from k — (tt, tt) by an amount 
larger than theoretically predicted. Additional confu- 
sion arose when experiments were performed in other 
materials. It was found that Nd2-xCexCu04 does not 
show incommensurability pO[ | while for Y Ba2Cu30Q+x 
(YBCO) the result is unclear. ||ll| Several authors sug- 
gested that the addition of a material dependent t' term 
to the Hubbard Hamiltonian could explain the quantita- 
tive, as well as the qualitative, differences observed. All 
previous studies of the U — t — t' model have been car- 
ried out in the weak coupling limit [U/t < 4). There are 
some disagreements among the published results. Benard 
et al. pjl^j concentrated on negative values of t' and used 
a perturbative approach. They found that incommensu- 
rability starts developing at a finite hole density which 
increases from half-filling as t' becomes more negative. 
A similar result, within RPA, was obtained in Ref. j|] 
though these authors expect incommensurability to dis- 
appear for t' = —0.45 due to a decreasing intensity of 
the peaks in the spectral function. Furukawa and Imada 
Ip^ studied positive values of t' using a zero temperature 
Monte Carlo algorithm for U/t = 4 and t' = 0.25. They 
found incommensurability, but with the peaks in S(k) 
moving along the diagonals in the Brillouin zone rather 
than towards the points X — {tt, 0) and Y = (0, tt) as 
found in the experiment and as it happens for t' = 0. 
Gooding et al. performed exact diagonalization stud- 
ies of the t — t' — J model that should resemble the strong 
coupling limit of the U — t — t' model. To increase the 
number of available momenta along particular directions 
in momentum space they did not use square clusters. For 
negative values of t' they observed incommensurability 
with S(k) moving along the diagonals, while for positive 
values of t' the peak in S(k) remained at (tt, tt) for all the 
dopings they studied. 

Due to the disparity of all these results and to the lack 
of numerical data on large square lattices at more realis- 
tic couplings larger than U/t = 4, we decided to carry out 
a systematic study of the structure factor at J7/t = 6 on 
8x8 lattices using quantum Monte Carlo. We worked at 
a temperature T/t — 1/4 since due to sign problems it is 
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not possible to study lower temperatures away from half- 
filling. Our results are presented in Fig. 2 where we show 
S(k) along the directions (0,0) - (tt, 0) - (tt, tt) - (0,0) 
{oi-r-X - M -T) for U/t = 6 at different fillings and 
for positive and negative values of t'. At half-filling, S(k) 
is independent of the sign of t' due to the symmetry in 
the Hamiltonian described in Sec. II. As the density is re- 
duced to (n) — 0.9, the strength of S(k) at the peak is 
slightly smaller for negative t' (Fig.2-a) than for positive 
t' (Fig.2-b). For negative (positive) hopping the maxi- 
mum in S(k) starts moving along the lines from (tt, tt) 
to (7r,0) and (0,7r) and it reaches k = (tt, 37r/4) and 
(37r/4,7r) for (n) = 0.76 (0.67). The figure clearly shows 
that the intensity of the peak in S(k) is drastically re- 
duced as the doping increases. In fact, when incommen- 
surability appears for positive t' (Fig.2-h) the peak can 
barely be seen. However, one has to remember that these 
features could be strongly enhanced as the temperature 
decreases. 

We conclude that negative values of t' seem to favor 
incommensurability more than positive values do. For 
both signs of t' we found that doping away from half- 
filling, S(k) peaks at k = (tt, fcmTr) and {kmi^^ir) where 
< fcm < 1, i-e. the peak does not move along the di- 
agonal contrary to the results of Ref. |^,|l3|. In Fig. 3 we 
plot km vs (n) for t' = —0.2, and 0.2 and, for compari- 
son, we also show the experimental results for LSCO 
It is clear that t' > 0(< 0) slightly decreases (enhances) 
incommensurability in comparison with the t' = case. 
The results for t' = —0.2 are closer to the experimental 
data but to reach a better agreement, it appears that 
larger negative values of t' should be used {t' w —0.4 
linearly extrapolating the results at t' — 0.2, and - 
0.2). This is an example where the value of t' suggested 
by comparing Fermi surface shapes obtained from band 
structure calculations with the ones obtained in the non- 
interacting case {t' = —0.2 for LSCO) does not provide 
a good fit for magnetic properties. 

We have tried to study the case t' = —0.45, which is 
the value suggested for YBCO. Due to sign problems, 
which become worse as \t'\ increases, we have been able 
to study only the case U/t = 4 and T = t/4 on an 8 x 8 
lattice. We found that the onset of incommensurability 
occurs at lower doping ((n) « 0.8) but the intensity of 
the peak in the structure factor, S(k), decreases and it is 
very difficult to decide whether incommensurability has 
developed. This resembles the experimental results for 
YBCO jll) that do not allow to determine if the peak 
in S(k) shifts away from (tt, tt). The authors of Ref. ||] 
found a similar behavior. ||lj] 

Previous analytical work |12|] predicted that km 
would become different from 1 immediately upon dop- 
ing for t' — but only when the density reaches a critical 
value (ric) < 1 for t' < 0. Since we work on an 8 x 8 lattice 
we only can detect incommensurability when km jumps 
from 1 to 0.75 and thus we cannot check the above men- 
tioned behavior. However, we have not observed quali- 
tative differences for the studied values of t'. Our results 



do not support the claims made in Ref. |^ where it was 
reported that a positive value of t' does not induce incom- 
mensurability. We found that, in this larger hole 
doping than was studied in Ref. (H is required to observe 
the effect but the tendency to incommensurability even- 
tually occurs. One of the reasons for the disagreement 
maybe the different symmetry of the clusters studied in 
Ref. [| that may induce incommensurability along the 
diagonals rather than towards X and Y as in the square 
clusters that we studied, and as it has been experimen- 
tally observed. 

To summarize, in this section we have studied the be- 
havior of the structure factor as a function of doping for 
positive and negative values of the diagonal hopping t' . 
We found that t' reduces the strength of the antiferro- 
magnetic correlations at half-filling. Away from (n) = 1 a 
negative hopping promotes incommensuration at a lower 
doping than a positive one. Using t' — —0.2 we did not 
find good agreement with experimental data for LSCO 
showing that deducing bare values of t' based on PES 
data may be misleading. For t' = —0.45 incommensura- 
bility was difficult to observe due to the strong reduction 
in the intensity of the peak in the structure factor. This 
behavior resembles experimental results for YBCO. 



IV. DENSITY VERSUS CHEMICAL POTENTIAL 

A very important property of the Hubbard model is the 
antiferromagnetic gap that appears at half-filling. This 
gap can be observed in the density of states by studying 
the spectral function A(k, oj) (see Sec. VII) or by analyz- 
ing the behavior of the density (n) vs fi, where fi is the 
chemical potential. It has been found that to change the 
doping from holes to electrons, ^ has to cross the gap 
and a plateau at (n) = 1 appears in the (n) vs n curve. 
[f7|,p^ Since the gap increases with U, so does the size of 
the plateau. In this Section we will study the influence 
of t' on these results. 

Before describing the numerical data let us discuss the 
expected behavior in the SDW mean field approximation. 
We found that for the range of parameters analyzed in 
this paper, the value of A that satisfies the mean field 
equations (see Appendix) changes very little when t' is 
introduced. For example, working at U/t— 4, A is 1.37 
for \t'\ = 0.2, i.e. only slightly smaller than at i' = 
where it takes the value 1.38. However, the main effect 
of t' is in a distortion of the bands in such a way that the 
actual gap in the density of states (DOS) is smaller than 
A. For the values of t' studied, this effect is important 
only when A is small. 

In Fig. 4 we show the mean field energy bands along 
F- M — A" — Fin momentum space, for t' —0 and —0.2 
and different values of A |l^. In Fig.4-a the bands for 
U/t = 6 are shown; the gap is equal to A which is 2.48 
and it is defined by the difference of energy between the 
two bands at k = S = (7r/2,7r/2) and k = X = {n^O) 
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(actually at t' = 0, the energies -E(k) are degenerate 
along the line from Y = (0, tt) to X). When t' = -0.2 
the effective gap has been reduced to 2.08 although A 
remains at 2.48, as can be seen in Fig.4-b, and it is then 
83% of the t' = value. Now the degeneracy along X — 
Y has been removed and the effective gap is defined by 
E{X) in the conduction band and -E(E) in the valence 
band. Results for U/t = 4 are presented in Fig.4-c and d. 
In this case the effective gap is reduced from 1.38 to 0.97, 
being now 70% of the original, when t' changes from to 
-0.2. It is clear that as U/t decreases the effect of t' on 
the size of the gap becomes more important. In Fig.4-e 
we present results for A = 0.5. This value is interesting 
because it is well known that the effective size of the 
gap is reduced in the weak-coupling case due to quantum 
fluctuations and this value of A could provide a more 
accurate representation of the U/t = 4 numerical data. 
Note that in this case, when t' = —0.2 the effective gap 
is very small. Then we would expect to see the plateau in 
{n) vs fi very much reduced. For U/t = 6, on the other 
hand, the gap changes only slightly at finite \t'\ — 0.2 
with respect to its value for t' = (see Fig.4-a and b). 
Then, as U /t increases, keeping t' fixed, our expectation 
is that the effect of t' on the size of the gap will become 
increasingly irrelevant. 

Now let us analyze the same problem but using QMC 
techniques. In Fig.5-a,b,c we present (n) vs /i for t' — 
and ±0.2 on an 8 x 8 lattice at f3t — 6 in the non- 
interacting case. When t' is negative (positive), half- 
filling is achieved for a negative (positive) value of the 
chemical potential rather than at /i = 0. It is clear from 
the figure that, as we pointed out in Section II, the curve 
for t' = 0.2 can be obtained from the curve at negative 
t' using the prescription (n)(t',/i) 2 — (n){—t',—^) 
even with interactions present. Then, below we will only 
present results for t' ~ —0.2. 

What is the effect of a finite Coulombic repulsion? In 
Fig.5-d and e we present (n) vs /i for U/t — 4 and (3t — Q 
on an 8 X 8 lattice. When t' = G (Fig.5-d), {n) has 
a plateau centered at 1 and it is clear that the chemi- 
cal potential has to move across the gap to change the 
doping from electrons to holes. For t' — —0.2 (Fig.5-e) 
the plateau is much smaller and the type of dopants can 
be changed from holes to electrons varying the chemical 
potential by a small amount. This behavior is in agree- 
ment with experimental data for the high Tc cuprates 
obtained by J. Allen's group [jl^ who observed that the 
chemical potential does not move across the gap as the 
dopants are changed from electrons to holes (although 
these results have not been yet reproduced using other 
experimental techniques). One alternative explanation 
of our numerical results may be that the ground state in 
this case is paramagnetic and the antiferromagnetic gap 
does not exist, but as we showed in Section III (Fig.l- 
a) long range antiferromagnetic order has already devel- 
oped in our finite cluster for the couplings we used thus 
this idea is incorrect. The SDW mean field calculation 
presented in the Appendix allows us to better interpret 



the data if we assume that the antiferromagnetic gap for 
U/t = 4 is smaller than the mean field prediction due 
to quantum fluctuations. If this is the case, the effective 
gap almost disappears as it was shown before in Fig.4-f. 
As the Coulombic repulsion is increased to U /t — Q the 
plateau at (n) = 1 reappears for t' = —0.2 as can be seen 
in Fig.5-g. In this case the plateau is similar to the one 
observed for t' = (Fig.5-f). This behavior is also in 
agreement with the SDW mean field results. The curves 
presented in Fig. 5 appear to be very smooth and they 
do not indicate the presence of phase separation in the 
system. 

Summarizing, in this section we studied the behavior 
of the density of electrons as a function of the chemical 
potential for positive and negative values of a diagonal 
hopping t' . We found that t' reduces the size of the gap. 
We also found that in weak coupling {U/t < 4) and at 
\t'\ — 0.2 the chemical potential does not move across 
the naive gap of order U when the dopants change from 
holes to electrons in agreement with the PES results of 
Allen et al. jl^ At larger U/t, the effect of t' becomes 
less important as observed in previous sections. 



V. LOCAL MAGNETIC MOMENTS AND 
DOUBLE OCCUPANCY 

The square of the magnetic moment per site is re- 
lated to the probability of double occupancy S — 
through the expression 

((m?)2) = (n)-2(niTnu). (3) 

In the non-interacting case the probability of dou- 
ble occupancy is independent of t' , and it is given by 
{n^riii) = (ri)^/4. The effect of the interaction U is to 
suppress the probability of double occupancy. In the 
interacting case {ni-^nj^) depends on t' . As it can be seen 
in Fig. 6, a negative t' decreases 5 when compared with 
t' ~ Q while a positive value increases it. At half- filling 
and using the symmetry described in Section II, it can 
be shown that the probability of double occupancy is 
independent of the sign of t' . According to Eq.(3) the 
mean-square magnetic moment per site will slightly in- 
crease for negative t' and decrease for positive t' . 

VI. PHOTOEMISSION EXPERIMENTS AND 
FERMI SURFACE 

In the last few years much progress has been made 
in the development of angle-resolved photoemission 
(ARPES) techniques and reliable data for the quasiparti- 
cle dispersion near optimal doping for Bi2Sr2CaCu20^ 
(Bi2212) and YBCO have been obtained. |]l| The study 
of the quasiparticle peak as a function of energy and 
momentum provides the shape of the Fermi surface for 
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the different materials. Until recently, some experiments 
seemed to provide evidence for large electron-like FS 
but these results have been lately challenged by Aebi et 
al. |2C[| Using a new photoemission technique these au- 
thors have reported the presence of a so-called shadow 
band probably due to antiferromagnetic correlations , 
as well as the existence of features resembling hole pock- 
ets [H in Bi2212 (see Fig. 7). ARPES also provides in- 
formation on the quasiparticle bandwidth which appears 
to be of the order of a fraction of eV independently of 
the material. Note that this bandwidth is much smaller 
than the one predicted by LDA calculations |l] which 
casts doubts on the practice of selecting values of t' with- 
out properly considering the strong correlations. Flat 
regions about X and Y have also been found in the hole- 
quasiparticle dispersion. Recently, Wells et al. ||2^ re- 
ported ARPES measurements on the antiferromagnetic 
insulator Sr2CuOCl2- This material is very difficult to 
dope but the behavior of the quasiparticle at half-filling 
is important as a probe for the microscopic theories for 
holes in antiferromagnets that have been proposed. The 
authors of Rcf. tried to fit their data with a t — J 
model with J = 0.125eV and found good agreement for 
the bandwidth but they could not fit the dispersion along 
the X-Y line. In a recent paper Nazarenko et al. |Q sug- 
gested that the addition to the t — J model of a t' hop- 
ping, which would be material dependent, could help to 
improve the fitting. Using the Born approximation they 
obtained agreement with experiments along the X-Y line 
introducing a.t' = — 0.35i, with J /t — 0.3 and J = 0.125 
eV. Along the X — T line disagreement with the exper- 
imental data remained although along this line a quasi- 
particle peak in the experiment is not clearly observed 
and the error bars are very large. Using the SDW mean 
field approximation described in the Appendix we calcu- 
lated the dispersion and we tried to fit the experimental 
data. Using U/t=l{) and t' = — 0.2t we found reasonable 
agreement with the experimental result as can be seen in 
Fig. 8. These results confirm the importance of t' terms 
to describe the cuprates. 

As it was stated in the Introduction, one of the rea- 
sons often invoked to introduce a t' hopping in models 
for high Tc superconductors is to match the shape of 
the Fermi surface to experimental results or to LDA cal- 
culations. Q The non-interacting Hubbard model has a 
nested Fermi surface (FS) at half-filling and it becomes 
large and electron-like immediately upon doping. The 
addition of a negative t' removes the nesting at half-filling 
and, instead, the FS becomes open and hole-like. Upon 
further doping it eventually closes but it retains hole-like 
characteristics. It becomes electron-like only when the 
hole doping has been further increased. ||] However, it is 
possible that the introduction of electronic interactions 
could modify substantially the non-interacting FS. For 
example, in the two dimensional Hubbard model, the 
uniform magnetic susceptibility always decreases upon 
doping in the non-interacting case, but it increases when 
interactions are strong enough (in agreement with ex- 



perimental results). Thus, it is very important to study 
how interactions affect the shape of the FS. 

Since we work on finite lattices it is difficult to ac- 
curately determine the FS of a model. Ideally, the 
FS should be obtained examining the spectral function 
yl(k, w). However, as we show in Section VII presently 
the spectral functions cannot be calculated with enough 
accuracy and alternative ways of determining the FS 
must be considered. In some previous papers the 
FS was obtained by calculating the position of the mo- 
menta where n(k) k, 0.5. This approach works well for 
weakly interacting systems, but it may cause problems 
when strong interactions are considered. One can imag- 
ine that for a system with a "hole pocket" centered at 
S = (7r/2,7r/2), n(k) along the (0,0) - (7r,7r) line will 
jump from a value larger than 0.5 to zero as the pocket is 
encountered increasing the momentum away from (0,0), 
defining a first Fermi surface, and then, at the second 
FS, it will increase again but to some value smaller than 
0.5. Then the ?T.(k) w 0.5 criterion is clearly incomplete 
in strong coupling because it may miss important struc- 
ture. We can avoid this problem by assuming that a 
FS is likely to exist where the numerically obtained ?T.(k) 
changes very rapidly. In fact, since experiments actually 
only detect crossovers in n(k) it is important to know 
where they occur in numerical simulations. Thus, a study 
of the values of n(k) will shed light on possible changes 
in the shape of the FS as shown below. 

Our mean field analysis predicts strong modifications 
in n(k) as U /t increases at a fixed t' . In particular, at 
half-filling, fixed t' /t and large U /t, n(k) resembles more 
the t' = Q case with the same U/t than the i' 7^ 0, 
U/t ~ case. This can be seen from Eq.(9) in the 
Appendix. If the gap between the conduction and the 
valence band is such that at half-filling the conduction 
band is completely filled and the valence band is empty, 
then ri(k) is given by the same analytical expression as 
for the t' = model. Indeed, this is what we have ob- 
served with QMC. In Fig.9-a we present n(k) along the 
line X = (tt, 0) to y = (0, tt), at half-fiUing for U/t = 0, 4 
and 6 and t' — —0.2. The mean field value along this line 
agrees with the exact result for t' = Q {n{k) — 0.5) for 
the finite values oi U/t studied. It appears that as U/t 
increases with respect to \t'\/t, n(k) becomes more and 
more similar to the t' = case, i.e. the concavity tends to 
disappear and the curve becomes flatter as U /t increases 
which is very different from the non-interacting curve as 
can be seen in the figure. This result at half-filling is not 
unexpected since as U/t increases, double occupancy de- 
creases, the diagonal hopping becomes less relevant and 
the spin part of the model becomes effectively described 
by a nearest neighbor only Heisenberg Hamiltonian (since 
J' = U'^/U = 0.04J « J for \t'\ = 0.2), i.e. the frus- 
trating J' term is negligible. 

Now let us study what happens away from half-filling. 
In Fig.9-b, c and d, we present n(\i) for the same param- 
eters as in Fig.9-a but for (n) = 0.9, 0.8 and 0.7. We also 
show the corresponding mean field results (open symbols) 
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for the interacting case. Mean field results are expected 
to become less reliable as (n) moves further away from 
1 and antiferromagnetic correlations decrease. However, 
the MC points are always closer to the mean field values 
than to the non-interacting ones. 

For negative t' , the low temperature mean field anal- 
ysis predicts the formation of hole pockets upon doping 
about the S = (7r/2,7r/2) point. On the other hand, 
in the non-interacting case the hole-like FS is centered 
about M = (tt, tt). In addition, in the mean field approx- 
imation at zero temperature there are crossovers in n(k) 
denoted by the crosses in Fig.lO-a where T = mean 
field results on a 100 x 100 lattice at (n) = 0.9, t' = -0.2 
and U/t = 6 are presented. The full squares define the 
actual FS and the crosses signal the points where n{k) 
changes the most rapidly. The doted line indicates the 
points in momentum space where n(k) = 0.5. Notice 
that Fig.lO-a resembles the results obtained by Aebi and 
collaborators, |20| (see Fig. 7) if we identify the points 
of maximum variation in n(k) with the experimental FS. 
As the temperature is increased within the mean field ap- 
proximation the Fermi surface around the pockets disap- 
pears and, instead, rapid crossovers are observed, i.e., 
the solid squares in Fig.lO-a are replaced by crosses. For 
U/t — 6 the pockets remain observable up to T/t « 1/8. 
As the temperature increases to T/t = 1/4, i.e., the value 
used in our Monte Carlo simulations, it can be seen in 
Fig.lO-b that the pockets about S have disappeared and 
only the lines of crossovers remain. The open squares in- 
dicate the points where n(k) — 0.5. Then, the n(k) = 0.5 
criterion applied to the MF data would incorrectly indi- 
cate a closed FS which does not resemble the actual FS 
at T = at all! In Fig.lO-c we show the FS at T = 
for the non-interacting system. In this case, there is a 
jump in n(k) from 1 to at the FS and the two crite- 
ria to determine its position agree. As the temperature 
increases to T = 1/4, Fig.lO-d, in the non-interacting 
case, both criteria still agree but a spurious crossover is 
induced along the direction r — X{Y) due to temperature 
effects. However, this variation of n(k) is much smaller 
than along the directions where a FS exists at T = 0. 

With this qualitative discussion in mind let us analyze 
our QMC results obtained on an 8 x 8 lattice at T — t/A 
using U/t = 6 and t' = —0.2. We will consider a sys- 
tem slightly doped away from half-filling with (n) = 0.9. 
The open squares in Fig.lO-e denote the points where 
n(k) « 0.5 while the crosses denote the regions of max- 
imum variation of n(k). These results are in excellent 
agreement with the mean field calculation on an 8 x 8 
lattice at the same temperature. Also notice that this 
figure resembles Fig.lO-b where mean field results at the 
same temperature on a 100 x 100 lattice are presented. 

Now let us point out the differences with the non- 
interacting case shown in Fig.lO-f which is equivalent 
to Fig.lO-e but with U=0. Notice that the criterion 
ri(k) « 0.5, denoted by the open squares, provides a 
"Fermi surface" identical to the one for the U/t = 6 case 
(Fig.lO-e) . However, the position of the maximum vari- 



ation of ".(k) (crosses) is remarkably different. [ p7[ 

Then, one of the main results of this paper is that for 
intermediate and large values of the Coulomb interaction 
it is incorrect to assume that the shape of the Fermi sur- 
face is similar to the non-interacting one. The criterion 
of maximum variation of n(k) does not support such an 
assumption. In this case, a better agreement is obtained 
when comparing numerical data with a SDW mean field 
calculation. According to this mean field approach, hole 
pockets about S should be observed at temperatures as 
high as T/t — 1/8. This result is in agreement with ex- 
periments for Bi2212 ||2^. Unfortunately, for U/t = 6 in 
the interesting doping regimes (0.8 < (n) < 0.9), we are 
not able to obtain accurate results at temperatures lower 
than T/t — 1/ A because the average value of the sign of 
the fermionic determinant becomes very small to directly 
verify the existence of hole pockets. 

As a final comment we want to add that a study of 
the points where n(k) varies more rapidly for the t' — 
Hubbard model provides a result qualitatively similar to 
the one displayed on Fig.lO-e suggesting that a negative 
t' term may be dynamically generated, For positive 
t' hole pockets appear at X, Y and symmetrical points. 

In this section we have measured n(k) numerically and 
we studied its behavior finding the points in momentum 
space where it changes more rapidly. The location of 
these points for negative values of t' determines a shape 
very similar to the one experimentally obtained in Ref. 
[EOl where the Fermi surface of Bi2212 was studied. 



VII. DYNAMICAL PROPERTIES 

It is very important to understand how the introduc- 
tion of t' affects the spectral function A(k, w). To eval- 
uate this quantity we used quantum Monte Carlo and 
the maximum entropy technique, In Fig. 11 we show 
A{k, Lu) for J7/i = 4 on an 8 X 8 lattice, at T = t/6, half- 
filling, and for t' = (Fig.ll-a) and -0.2 (Fig.ll-b). At 
half filling and t' = the chemical potential is located 
inside the gap. A(k, lo) is the same along the line from 
(0, tt) to (tt, 0) in agreement with the SDW mean field 
result and with the numerical observation showing that 
ri(k) — 0.5 along this line. The size of the gap agrees 
with the result displayed in Fig.5-d. When t' = —0.20, 
as discussed in Section IV, the gap is reduced and, within 
the accuracy of the maximum entropy procedure, only a 
pseudogap is observed (Fig.ll-b). This again, is in agree- 
ment with the results presented in Fig.5-e and with the 
SDW mean field calculations. Also notice that the chem- 
ical potential does not lie inside the pseudogap. It is clear 
that the symmetry along X — Y has been removed. 

The study of the spectral function should, ideally, de- 
termine the shape of the Fermi surface. This is done 
by observing the momenta where the quasiparticle peak 
crosses the chemical potential. The problem that we have 
observed with this approach is that the quasiparticle peak 
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has a finite width which is rather large due to lack of ac- 
curacy of the maximum entropy technique. In fact, the 
value of k at which the peak starts moving across the 
chemical potential from the left is very different for the 
value of k when the peak has finally moved completely to 
the right of the chemical potential, and this introduces 
large error bars. In Fig.ll-c we present the spectral func- 
tion for U/t = 8, t'/t = -0.2 at (n) = 0.9 and T/t =1/2 
on an 8 X 8 lattice. Along the direction F — X{Y) it 
can be seen that a broad quasiparticle peak moves across 
the chemical potential; this would suggest the existence 
of a Fermi surface in this direction as Aebi et al. |2^] 
found but we cannot determine accurately at what value 
of the momentum the actual crossing occurs. Along the 
line X{Y) — M there are indications that a Fermi sur- 
face may exist very close to X{Y). These Fermi surfaces 
would appear in the regions where we detected the largest 
variations of n(k) (see Fig.lO-e). The determination of 
the FS through the study of the spectral function was 
used in Ref. where the Hubbard model with U/t = 8 
and (n) = 0.87 was studied at T/t = 1/2 on up to 12 x 12 
lattices. The conclusion is that, in this case, there is a 
closed Fermi surface centered at k = (tt, tt). In Fig. 12 we 
present the shape of the Fermi surface that we obtained 
studying the problem on an 8 x 8 lattice and using the 
criteria described in Section VI. The n(k) w 0.5 criterion 
(open squares) gives an electron-like Fermi surface closed 
about k = (0, 0) as previously observed in Ref. but we 
have provided enough evidence that this criterion is not 
necessarily reliable. The points of rapid crossover in rt(k) 
(stars) provide two surfaces, similar to the behavior pre- 
dicted by the SDW mean field (see Appendix and Section 
VI) and the experimental data of Ref. |Q . 

In the figure we also plot the results obtained from 
Fig.l in Ref. [^. The full circles indicate the values of 
k where the maximum of the quasi-particle peak crosses 
the chemical potential and the error bars correspond to 
the width of the peak. Once the data of Ref. are 
supplemented by the proper error bars, the results are 
consistent with those obtained with the maximum vari- 
ation of n(k) criterion. However, due to the large error 
bars it is still not possible to decide whether the Fermi 
surface is closed about F or M. 

In this section we studied how a diagonal hopping t' 
added to the Hubbard Hamiltonian affects the behav- 
ior of the spectral function A(k, w). We observed that 
in weak coupling [U/t—^) the gap that appears at half- 
filling is reduced in agreement with the results presented 
in Section IV. We found that the data are not accurate 
enough to allow the unique determination of the shape 
of the Fermi surface. However, in the strong coupling 
regime {U/t = 8) away from half- filling, we found indica- 
tions of a Fermi surface consistent with the one presented 
in Section VI. 



VIII. CONCLUSIONS 

In this paper we analyzed the 2D Hubbard model with 
a diagonal hopping term in the intermediate coupling 
regime using Monte Carlo techniques. We observed that 
a diagonal hopping t' decreases the strength of spin-spin 
correlations. A negative diagonal hopping accelerates the 
onset of incommensurability while a positive one retards 
it. The intensity of the peak in the structure factor de- 
creases rapidly as the system is doped away from half- 
filling but its strength may increase at lower tempera- 
tures. Comparing our numerical results with experimen- 
tal data for LSCO, we found that a t' more negative than 
the suggested value |l] -0.2 appears to be needed to reach 
quantitative agreement. 

We found that t' tends to reduce the size of the anti- 
ferromagnetic gap but this effect is important only in 
the weak coupling regime. Double occupancy is en- 
hanced (reduced) by a positive (negative) diagonal hop- 
ping. This effect does not occur in the non-interacting 
case. 

A study of n(k) allows us to define a new criterion to 
find probable locations of Fermi surfaces working numer- 
ically at finite temperature. Since experimentalists also 
work at finite temperature and find the FS as the points 
in momentum space where n(\<i) has rapid crossovers, we 
looked at these points in our numerical data. We ob- 
served that lines of rapid crossovers in n(k) appear in 
positions similar to the FS obtained experimentally at 
room temperature for Bi2212 An analogous result 
is obtained for the t' — Q case which indicates that a 
negative hopping term may be dynamically generated. 

A study of the spectral functions show how the addi- 
tion of a diagonal hopping term removes the symmetry 
along the line X — Y and reduces the size of the gap at 
half-filling. 

We conclude that it is very dangerous to extract a bare 
parameter of the Hamiltonian like t' from experimental 
(PES) data where renormalizaed parameters play the im- 
portant role. 
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X. APPENDIX: SDW MEAN FIELD 



The spin density wave (SDW) mean field formalism has 
been applied to the Hubbard model and also to the 
U — t — t' model [^,0 . At half-filling the ground state for 
the U — t — t' model is antiferromagnetic it U > Uc, the 
critical coupling Uc is a function of t' and for \t'\ = 0.2 we 
found numerically Uc ~ 2.5. Since we will study U=4 
or higher we want to concentrate in the antiferromagnetic 
solution. In this case we found two energy bands given 

by 



— —At'coskxCosky — /i ± 



where 



A2, 



Efc = —2t{coskx + cosky), 



(4) 



(5) 



(6) 



and A is obtained by solving the self-consistent equa- 
tions: 



U ^ 



El 



(7) 



and 



(8) 

k 

where j{x) is the fcrmi function given by and 
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Figure Captions 

1. (a) Spin-spin correlation C(r) = (S?S?+r>(-l)l''l for 
U/t = 4, T t/6 on an 8 X 8 lattice at half-filling 
for <' = (open squares) and t' = —0.2 (filled 
squares). The error bars are of the size of the dots; 
(b) Same as (a) for U/t = 6, T = t/4. 

2. (a) Structure factor S(k) along the lines V — X — 
M - r for U/t = 6 and /3t = 4 on an 8 x 8 lattice 
at (n) = 0.9 and t' — —0.2; (b) Same as (a) for 
t' = 0.2; (c) Same as (a) for {n) = 0.8; (d) Same 
as (c) for t' — 0.2; (e) Same as (a) for {n) — 0.7; 

(f) Same as (e) for t' — 0.2; (g) Same as (a) for 
{n) = 0.6; (h) Same as (g) for t' = 0.2. 

3. Value of k„i as a function of doping for t' — 
0.2 (open squares), (triangles) and -0.2 (dark 
squares) for J7/t = 6 on an 8 x 8 lattice and j3t — 4. 
The crosses are experimental data from Rcf. 

4. (a) SDW mean field energy bands on an 8 x 8 lattice 
along V — M — X — T vn. momentum space, for A — 
2.48 that corresponds to U/t — 6 for t' ~ 0; the 
effective gap is shown with dashed lines; (b) Same 
as (a) but for t' — —0.2; (c) Same as (a) but for 
A = 1.38 that corresponds to U/t = 4. (d) Same 
as (c) but for t' = —0.2; (e) Same as (a) but for 
A = 0.5. (f) Same as (e) but for t' = -0.2. 

5. (a) The density (n) as a function of /i on an 8 x 8 
lattice for U/t = 0, T = t/& and t' = 0. (b) Same 
as (a) for t' = -0.2; (c) Same as (a) for t' = 0.2; 
(d) Same as (a) for U/t = A; (e) Same as (b) for 
U/t = 4; (f) Same as (a) for U/t = 6 with T = t/4:; 

(g) Same as (f) for t' = —0.2. The error bars are 
smaller than the symbols. 

6. Probability of double occupancy S as a function of 
filling (n) on an 8x8 lattice at /? = 4 and t = 1 for 
U=6 and t' = 0.2 (open squares), (crosses) and 
-0.2 (filled squares). The filled triangles are results 
for U/t==0. 

7. Fermi surface for Bi2212 as shown in Ref. | pO[ |. 



8. Quasiparticle dispersion of the U — t~t' model cal- 
culated using a SDW mean field using t' — —0.2t, 
U/t = 10 and U 4eV (sohd hne), t - J resuhs 
(dotted line) and experimental ARPES results from 
Ref. (21 . Through the relationship J = U'^/U 
we found J/t ^ 0.4 and using t — 0.4eV we get 
J = 0.16eV. 

9. Quantum MC values for n(k) along X — F on an 
8x8 lattice at /3t = 4 and t' = -0.2 for [/ = 
(crosses), U = 4 (filled circles) and U — 6 (filled 
squares) and mean field results indicated by open 
circles {U = 4) and open squares {U = 6) and 
connected with dashed lines for a) (n) = 1.0 (in this 
case the MF results along the direction shown are 
independent of U and open circles only are used), 
b) (n) = 0.9, c) {n) = 0.8, and d) (n) = 0.7. 

10. a) Mean field results for U/t = 6,T = and (n) = 
0.9 on a 100 x 100 lattice. The dark squares denote 
the Fermi surface; the crosses indicate the values of 
k where n(k) changes most rapidly and the dotted 
line indicates the values of k where ri(k) — 0.5; b) 
same as a) for T/t = 1/4. Here, there is no real 
Fermi surface and the values of k where ?T.(k) = 0.5 
are indicated by open squares; c) same as a) for 
U/t = 0;d) same as b) for U/t = 0; e) Monte Carlo 
resuhs for U/t = 6, T/t =1/4 and (n) = 0.9 on 
a 8 x 8 lattice. The crosses indicate the values of 
k where n(k) changes most rapidly and the open 
squares indicate the values of k where n(k) ~ 0.5;f) 
same as e) for U/t = 0. 

11. a) A{k,uj) for U/t = A,t'=^ 0, T/t = 1/6 on an 8x8 
lattice at half-filling. The moments are in units of 
7r/4. The dotted line indicates the position of the 
chemical potential; b) same as a) for t' = —0.2; c) 
A{k,uj) for U/t = 8, t' ^ -0.2, T/t 1/2 on an 
8x8 lattice at (n) = 0.9. 

12. Fermi surface for the 2D Hubbard model with 
U/t = 8, (n) = 0.87 at T/t = 1/2 using the 
n(k) — 0.5 criterion (open squares and solid line), 
the maximum variation of n(k) criterion (stars and 
dashed line) and from the maximum entropy re- 
sults of Ref. [^ (solid circles, with error bars, and 
dotted line). 
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Max CO = 10.0 

Max A(co) = 0.4 



